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Abstract. A simple model of oscillator chain with dynamical traps and additive white noise is considered. 
Its dynamics was studied numerically. As demonstrated, when the trap effect is pronounced nonequilibrium 
phase transitions of a new type arise. Locally they manifest themselves via distortion of the particle 
arrangement symmetry. Depending on the system parameters the particle arrangement is characterized by 
the corresponding distributions taking either a bimodal form, or twoscale one, or unimodal onescale form 
which, however, deviates substantially from the Gaussian distribution. The individual particle velocities 
exhibit also a number of anomalies, in particular, their distribution can be extremely wide or take a quasi- 
cusp form. A large number of different cooperative structures and superstructures made of these formations 
are found in the visualized time patterns. Their evolution is, in some sense, independent of the individual 
particle dynamics, enabling us to regard them as dynamical phases. 



PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.45.- 
linear dynamics and nonlinear dynamical systems - 05.70.Fh Phase transitions: general studies 
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1 Introduction 



For the last several decades various phenomena cause by 
the ordering action of noise in nonequilibrium systems are 
found (for a general review see Refs. m^lEE])- Popular 
examples are stochastic resonance |51|H1I7| , coherence reso- 
nance I21IH1, noise-induced transport and noise- induced 
phase transitions J. 10, 11^. Typically the latter are due to 
multiplicative noise. However, additive noise in the pres- 
ence of other multiplicative noise can also induce phase 
transitions |12|[T3.. 14 or individually cause a hidden phase 
transition to become visible jl5| . 

The constructive role of noise is peculiar to nonequi- 
librium systems only. In thermodynamic systems, for ex- 
ample, the phase formation is solely due to a certain reg- 
ular "force" changing its form, in particular, the number 
of stationary points. Available noise mainly perturbs the 
system motion around these points. 

In the present paper we pay attention to a new class 
of nonequilibrium systems, namely, many particle ensem- 
bles where some long-lived cooperative states can form 
whereas the regular component of "individual" forces has 
no stationary points except for one corresponding to the 
homogeneous state. The latter, in addition, is locally sta- 
ble for all the values of the system parameters. Besides, 
only additive noise enters such systems. 



Originally investigation of the model under considera- 
tion was stimulated by a wide class of intricate coopera- 
tive phenomena found in the dynamics of vehicle ensem- 
bles moving on highways, motion of fish and bird swarms, 
stock markets, etc (for a review see Ref. ^j). The model 
background is the following. People as elements of a cer- 
tain system cannot individually control all the governing 
parameters. Therefore one chooses a few crucial parame- 
ters and focuses on them the main attention. When the 
equilibrium with respect to these crucial parameters is at- 
tained the human activity slows down retarding, in turn, 
the system dynamics as a whole. For example, in driv- 
ing a car the control over the relative velocity is of prime 
importance in comparison with the correction of the head- 
way distance. So under normal conditions a driver, first, 
should eliminate the relative velocity between his car and 
a car ahead and only then optimize the headway. 

These speculations have led us to the concept of dy- 
namical traps, a certain "low" dimensional region in the 
phase space where the main kinetic coefficients specify- 
ing the time scales of the system dynamics become suffi- 
ciently large in comparison with their values outside the 
trap region |17II18II19II20| . As a result long-lived states 
have to appear. In time patterns these states manifest 
themselves like a sequence of fragments within which at 
least one of the phase variables remains approximately 
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constant. These fragments are continuously connected by 
sharp jumps of the given variable. Paper |18) demonstrated 
that such long-lived states do exist in the dense traffic 
flow and proposed some model of dynamical traps to ex- 
plain the observed features of car velocity time series. Pa- 
pers [T^ll^ simplified this model to single out the dy- 
namical trap effect in its own record. Similar phenomena 
seem to be observed in physical systems also, for example, 
during the non-monotonic relaxation of Pd-metal alloys 
charged with hydrogen [23 • Paper studied a single 
oscillator with dynamical traps and demonstrated numer- 
ically that white noise can cause the distribution function 
of oscillator position to convert from the unimodal form to 
the bimodal one. It is due to the fact that inside the trap 
region the regular "force" is depressed only rather than 
changes the sign and the system motion is mainly caused 
by a random Langevin "force" . A first step towards of this 
effect in oscillator ensembles was made in Ref . |201 ■ In par- 
ticular, the dynamical traps were demonstrated to be able 
to give rise to the system instability and an anomalous 
velocity distribution like a cusp oc exp{ — smoothed, 
naturally, inside a narrow transition region. It should be 
noted that similar anomalous velocity distributions were 
found for dense traffic flow j22j . 

The purpose of the present paper is to demonstrate 
that an ensemble of such oscillators with dynamical traps 
can exhibit a number of anomalous cooperative phenom- 
ena, their detailed investigation will be published else- 
where. However, first, we clarify the relation between the 
specific mathematical form of the model to be studied and 
the concept of dynamical traps discussed above. 
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Fig. 1. The particle ensemble under consideration and the 
structure of the phase space. The darkened region depicts the 
points where the dynamical trap effect is pronounced. For the 
relationship between the variables Xi, Vi, hi, and "di see formu- 
lae (21 and 0. 



retard the correction of its relative position. This assump- 
tion leads to the appearance of some common cofactor 
J7('i9j) in the governing equation like this 



oc ~f2{iSi)(rit + a-di 



The cofactor i7(i?) is to meet the inequality ^{1)) <C 1 
for -d -de and « 1 when -d 3> i^c, where -de is a 

certain critical value quantifying the particle perception 
of speed. Exactly the appearance of such a factor is the 
implementation of the dynamical trap effect. Now let us 
specify the model. 



Motivated behavior of particles 

Keeping in mind the aforesaid about the human behavior 
we consider a one-dimensional ensemble of "lazy" particles 
characterized by their positions and velocities {xi,Vi} as 
well as possessing some motives for active behavior. Par- 
ticle i wants to get the "optimal" middle position between 
the nearest neighbors. So one of the stimuli for it to accel- 
erate or decelerate is the difference rji — Xi~'^{xi^i+Xi+i) 
provided its relative velocity di = Vi — ^(vi-^i + Vi^i) with 
respect to the pair of the nearest neighbors is sufficiently 
low. Otherwise, especially if particle i is currently located 
near the optimal position, it has to eliminate the relative 
velocity, being the other stimulus for particle i to change 
its state of motion. Since a particle cannot predict the dy- 
namics of its neighbors it has to regard them as moving 
uniformly with the current velocities. So both the stim- 
uli are to determine directly its acceleration dvi/dt. The 
model to be formulated in the next section combines both 
of these stimuli within a linear approximation similar to 
{rji + (J'di), where a is the relative weight of the second 
stimulus. 

When, however, the relative velocity di of particle i 
attains sufficiently low values the current situation for it 
cannot become worse, at least, rather fast. So in this case 
particle i prefers not to change the state of motion and to 



2 Model 

The following linear chain of N point-like particles is con- 
sidered (Fig.^. Each internal particle i 1,N can freely 
move along the a:-axis interacting with the nearest neigh- 
bors, namely, particles i — 1 and i + I via ideal elastic 
springs with some quasi- viscous friction. The dynamics of 
this particle ensemble is governed by the collection of cou- 
pled equations 



— ^ = -f2{d,, h,)[r], + ad, + a^v,] + ei,{t) . (2) 
at 

Here for z — 2, 3, . . . , — 1 the variables r]i and "di to 
be called the symmetry distortion and the distortion rate, 
respectively, are specified as 

r]i= ^i- ^(^i-i + ' (^) 
i^i = Vi - ^(^i-i + ^i+l) ' (4) 

the mean distance hi between the particles at the point 
Xi, by definition, is 

hi ^ l-{xi+i - Xi-i) , (5) 
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and {^i(t)} is the collection of mutually independent white 
noise sources of unit amplitude, i.e. 



(6) 



In addition, the parameter e is the noise amplitude, a is 
the viscous friction coefficient of the springs, ctq is a small 
parameter that can be treated as a certain viscous friction 
related to the particle motion with respect to the given 
physical frame. It is introduced to prevent the system mo- 
tion as a whole with infinitely high velocity. Besides, the 
symbol (...) denotes averaging over all the noise realiza- 
tions, Sii' and 5{t — t') are the Kronecker symbol and the 
Dirac 5- function. The factor fl{di,hi) is due to the effect 
of dynamical traps and actually following our previous pa- 
per J7] the Ansatz 



l?2 + 1 



(7) 



with a function A(/i) such that 



boundary conditions © are fulfilled by assuming the wave 
number k to take the values km — ■nm/[{N — for 
m = ±1, ±2, . . . , ±(A^ — 2). For large values of the particle 
number A'^ the parameter k can be treated as a continu- 
ous variable. Using the standard technique the system of 
equations (0) for perturbation ifTTl) leads us to the fol- 
lowing relation of the instability increment 7 (A:) and the 
wave number k: 



7 = -■^o 



1 



-do + cr sm ^2 



kl 



'2Q0 sir? (y) - [21 i(7o + o- sin^ 



(12) 



In deriving expression 1121) Ansatz Q has been used, en- 
abling us to set J7o — [2{0,l) — Whence it follows 
that Re7(fc) > for fc > 0, so the homogeneous state of 
the chain is stable with respect to infinitely small pertur- 
bations of the particle arrangement. 



A\h) = + (1 - A 



ft2 
"0 



/l2 + 



hi 



(8) 



is used. The parameter A G [0.1] quantifies the dynam- 
ical trap influence and the spatial scale ho specifies the 
small distances within which the trap effect is to be de- 
pressed, i.e. for h <^ ho the value A{h) « 1 whereas when 
h ^ ho/ A the value A{h) « A. If the parameter A = 1, 
the dynamical traps do not exist at all, in the opposite 
case, A <C 1, their influence is pronounced inside a certain 
neighborhood of the h-axis (trap region) whose thickness 
is about unity (Fig. The temporal and spatial scales 
have been chosen so that the thickness of the trap region 
be about unity as well as the oscillation circular frequency 
be also equal to unity outside the trap region. The termi- 
nal particles, i — 1 and i = N, are assumed to be fixed, 
i.e. 

xi(i)=0, XNit) = iN-l)l, (9) 

where / is the particle spacing in the homogeneous chain. 
The particles are treated as mutually impermeable ones. 
So when the coordinate Xi and Xi+i of an internal particle 
pair become identical the absolutely elastic collision is as- 
sumed to happen, i.e. if Xi{t) = a;i+i(t) at a certain time 
t then the timeless velocity exchange 



v^{t + 0) = v,+i{t~0), 
+iit + 0) = v,{t - 0) 



(10) 



comes into being. The multiparticle collisions are ignored. 

The system of equations ((T|l- (|10(l forms the model un- 
der consideration. 

The stationary point xf = {i — 1)1 is stable with re- 
spect to small perturbations. It steams from the linear sta- 
bility analysis with respect to perturbations of the form 

6xi{t) (X exp{jt + ikl{i - 1)} , (11) 

where 7 is the instability increment, k is the wave num- 
ber, and the symbol i denotes the imaginary unite. The 



3 Nonlinear dynamics 

The nonlinear dynamics of the given system has been an- 
alyzed numerically. The integration of the stochastic dif- 
ferential equations (O was performed using the E2 
high order stochastic Runge-Kutta method |^ (see also 
work |24)). Particle collisions were implemented analyz- 
ing a linear approximation of the system dynamics within 
one elementary step of the numerical procedure and find- 
ing the time at which a collision has happened. Then this 
step treated as a complex one was repeated. The inte- 
gration time step of 0.02 was used, the obtained results 
had been checked to be actually stable with respect to de- 
creasing the integration time step. The ensemble of 1000 
particles was studied in order to make the statistics suf- 
ficient and to avoid a strong effect of the boundary con- 
ditions. The integration time T was chosen from 5000 to 
8000 time units in order to make calculated distributions 
stable. At the initial stage all the particle were distributed 
uniformly in space whereas their velocities were randomly 
and uniformly distributed within the unit interval. 

The results of numerical simulation were used to eval- 
uate the following partial distributions 



1 



N-M 



[N ~2M){T -To) 



J2 j dtS{z-z,{t)), (13) 



where the time dependence z^ (t) describes the dynamics of 
one of the variables r]i(t), '&i{t), and Vi(t) ascribed to par- 
ticle i and z is a given point of the space Rz describing the 
symmetry distortion 77, the distortion rate and the par- 
ticle velocity v, respectively. The variables {77, -d, v} enable 
one to represent the system dynamics portrait within the 
space M,, x x R„ or its subspace. Besides, N is the total 
number of particles in the ensemble and M is the num- 
ber of particles located near each of its boundaries. They 
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Fig. 2. The distribution functions of the coordinate rj and the velocity 1} of the movable particle in the 3-particle ensemble. 
These distributions were obtained by averaging the simulation results over time interval of 100000. Solid lines correspond to the 
case of strong trap effect, A = 0.1, dotted line match the absence of dynamical traps, A = 1.0. The right widows depict some 
path fragment formed by the movable particle during time interval of IQQQ time units. Other used parameters are e = 0.1 and 
two values of the dissipation rate a + ao ~ 0.1 (upper row) and a + ao = 1.0 (lower row). 



are excluded from the consideration in order to weaken 
a possible effect of tlie specific boundary conditions. The 
same concerns the lower boundary of time integration Tg, 
its value is chosen to eliminate the effect of the specific 
initial conditions. 

The numerical implementation of the integration over 
time in expression was related to the direct summa- 
tion of the obtained time series and the partition of the 
corresponding space was chosen so that the results be 
practically independent of the cell size. The value of M 
was also chosen using the result stability with respect to 
the double increase in M. Typically the value M ~ 50 was 
chosen for N = 1000, for = 3, naturally, M = 1, and 
To - 500-1000. 



3.1 Three particle ensemble 

The given oscillator chain made of three particles is ac- 
tually the system studied in part previously jl7| . In this 
case only the middle particle is movable and the variables 
V '■= ^2 and "d := '62 are its coordinate and velocity. Here 
we also present the results for the 3-particle ensemble in 
order to have a feasibility of distinguishing characteristics 
of local nature from many particle effects. 

Figure[21compares the distribution functions V{rf) and 
V^d) obtained in the cases where the dynamical trap effect 
is absent (A — 1) and when the dynamical traps affect 
the particle motion substantially (A ^ 1). The upper 
windows correspond to the system with weak dissipation. 



(7 + (To = 0.1, whereas the lower ones are related to the 
case of strong dissipation, cr + ctq = 1.0 

In agreement with the previous results ^\ it is seen 
that the decrease of the parameter A, i.e. the dynamical 
trap intensification induces the conversion of the function 
7^(77) from the unimodal form to the bimodal one, with 
the dissipation no more then weakening this effect. A new 
result is the essential dependence of the velocity distribu- 
tion on the dissipation rate. In the case of weak dissipation 
the movable particle performs alternatively fast motions 
outside the trap region and slow motion inside it. The 
fast motion paths connect the neighborhoods Q_, Q+ of 
the 'P(77)-function maxima, whereas the slow motion arises 
when the particle wanders inside these regions. This fea- 
ture is visualized in Fig. |21 the right upper window shows 
a fragment of the particle path of duration about 1000 
time units. Therefore the obtained distribution function 
7^(1?) as seen in Fig. |21 (middle upper window) actually 
is made of two monoscale components. For the case of 
strong dissipation the two neighborhoods Q_ and (Q+ are 
not directly connected by the fast motion paths (Fig. |21 
right lower window). Now they rather uniformly spread 
over a certain domain on the {"q, z?}-plane, previously, they 
were located inside a sufficiently narrow layer. As a result 
the velocity distribution converts into a monoscale func- 
tion having a quasi-cusp form oc exp{ — It?!}. We relate the 
cusp formation to the properties of the system dynam- 
ics near the trap region. It is justified in Fig. 13 showing 
the resulting velocity distribution of the movable particle 
when noise and dissipation are absent. In this case the 3- 
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-1.0 -0.5 0.0 0.5 1.0 

middle particle velocity, u 



Fig. 3. An example of the velocity distribution 'P{'&) formed 
by the movable particle of the 3-particle ensemble without 
noise and dissipation (e = and cr = 0). The path on the 
phase plane {rj, i}} formed by this particle is show in the in- 
cluded window. In numerical simulation A = 0.1 was used and 
a phase path of the velocity amplitude about unity was chosen. 



particle system admits conservation of a certain "energy" 
and the phase paths form a collection of closed curves on 
the {77, z?}-plane Q^- 

If the dynamical trap effect is absent, A = 1, all these 
distributions, as it must, are of the Gaussian form shown 
in Fig. 13 with dotted lines. 

3.2 Multi-particle ensemble 

To analyze cooperative phenomena arising in such systems 
the dynamics of 1000-particle ensembles was implemented. 
Let us, first, consider local properties exhibited by these 
ensembles. The term "local" means that the correspond- 
ing state variable can take practically independent values 
when the particle index i changes by one or two. The vari- 
able Tji (expression 0) may be regarded in such a man- 
ner. It describes the symmetry of particle arrangement in 
space, when rji = particle i takes the middle position 
between the nearest neighbors, particles i — 1 and i -|- 1. 
A nonzero value of 77^ denotes its deviation from this po- 
sition, in other words, a local distortion of the ensemble 
symmetry. The latter was the reason for the used name of 
the variables rji as well as the variables iJ^ — drji/dt. 

Figure 0] exhibits the found distributions of the vari- 
ables rj and d depending on the dissipation rate a and 
the initial distance I between the particles, i.e. their mean 
density. Comparison of Fig. Eland Fig. 01 shows us that in 
this case of weak dissipation the distribution functions of 
the symmetry distortion 7^(77) and the distortion rate 
are qualitatively similar to those of the corresponding 3- 
particle ensemble. Only a few new features appear. First, 
for the system with high particle density (Z = 5) a small 
spike is visible at the center, 77 = 0, of the distribution 
function ^(77), which is pronounced in the case of strong 
dissipation. It corresponds to the symmetrical state of the 
particle ensemble being stable without dynamical traps 
and destroyed for the 3-particle ensemble. In the given 



case "many-particle" effects seem to reconstruct it in part. 
So in the given case the particle arrangement is charac- 
terized by three states, two of them match the extrema of 
the distribution function 'P{ri) and the symmetrical state 
singled out to a some degree. 

As for the 3-particle ensemble the distortion rate dis- 
tribution is again composed of two monoscale components, 
narrow and wide ones. Previously we have related them 
to the fast and slow motions. Figure 21(upper second win- 
dow) also justifies this. The narrow component is due to 
the particle motion inside the trap region and should be 
practically independent of the mean distance between par- 
ticles. By contrast, the wide one is to depend remarkably 
on the particle density because it matches the fast motion 
of particles outside the trap region and, thus, has to be 
affected by their relative dynamics. Exactly this effect is 
demonstrated in Fig. 01 visualizing also the corresponding 
properties of the particle paths. 

For the 1000-particle ensemble with strong dissipa- 
tion, CT w 1.0, the situation changes dramatically, although 
the characteristic scales of the corresponding distributions 
turn out to be of the same order in magnitude. In the 
given case the distribution function Virf) of the symmetry 
distortion has only one maximum at 77 = 0, however, its 
form is to be characterized by two scales. In other words, 
it looks like a sum of two monoscale components. One of 
them is sufficiently wide, its thickness is about the same 
value that is obtained for the corresponding particle en- 
semble with weak dissipation. Exactly this component ex- 
hibits a remarkable dependence of the particle density, 
enabling us to relate it to the particle motion outside the 
trap region. The other is characterized by an extremely 
narrow and sharp form shown in detail in the inner win- 
dow in Fig.01for the dense particle ensemble. Its sharpness 
leads us to the assumption that "many-particle" effects in 
such systems with dynamical traps cause the symmetrical 
state to be singled out from the other possible states in 
properties. 

By contrast, the distortion rate behaves rather sim- 
ilar to the previous case except for some details. When 
the mean particle density is high {I = 5) the wide com- 
ponent of the distortion rate distribution disappears and 
only the narrow one remains, with the latter having a 
quasi-cusp form oc exp{ — 17?|}. For the system with low 
density the peak of the distortion rate distribution splits 
into two small spikes. 

These features can be explained by applying to the low 
row right windows in Fig. ^ which exhibit typical path 
fragments formed by motion of a single particle on the 
{77'i?}-plane. Roughly speaking, now three motion types 
can be singled out: some stagnation inside a narrow neigh- 
borhood of the origin {77 = , i9 = 0} (visible well in the 
right window), slow wandering inside the trap region that, 
on the average, follows a line with a finite positive slope 
(clearly visible in the left window), and the fast motion 
outside the trap region (visible again in the left window). 
The fast motion fragments typically stem from an arbi- 
trary point of the low motion region and lead to a certain 
neighborhood of the origin. It seems that for the system 
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Fig. 4. The distribution functions of the symmetry distortion rj and the distortion rate i9 for the 1000 particle ensemble with 
low {I = 50, label 1) and high {I = 5, label 2) density and weak (a ~ 0.1) and strong (a ~ 1.0) dissipation. The right four 
windows depict characteristic path fragments of duration of 1000 time units formed by a single particle with index i = 500 
on the phase plane {rj, i}} which was chosen due to its middle in the given ensemble. The other used parameters are the noise 
amplitude e — 0.1, the trap effect measure A = 0.1, the small regularization friction coefficient ao = 0.01 and the regularization 
spatial scale ho = 0.25. The time interval within which the data were averaged changed from 2000 up to 5000 in order to make 
the obtained distributions stable. 



with low density particles have possibility to go sufficiently 
far from the origin and during the fast motion come into 
the stagnation region rarely. As a results, first, the dis- 
tortion rate distribution function is of a two scale form 
and contains two spikes on the peak. In the case of high 
density the fast motion is depressed substantially and the 
system migrates mainly in the slow motion region enter- 
ing the stagnation region for many times. So the distortion 
rate distribution converts into a single-scale function and 
the appearance of the symmetric state becomes often, giv- 
ing rise to a significant sharp component of the distortion 
distribution located near the point rj = 0. 

Now we discuss the nonlocal characteristics of the 100- 
particle ensembles. Figure |S1 depicts the velocity distribu- 
tions. As is seen it depends essentially on both the param- 
eters, the mean particle density and the dissipation rate. 
When the mean particle density is low and the dissipation 
is weak {I = 50 and a w 0.1) the velocity distribution is 
practically of the Gaussian form, however, its width gets 
extremely large values about 10. We recall that without 
dynamical traps the width of the corresponding distribu- 
tion does not exceed 0.5 (Fig. The tenfold increase of 
the particle density, I : 50 <—> 5,, shrinks the velocity distri- 
bution to the same order and its scale get values similar 
to that of the distortion rate distribution in magnitude. 
However in this case the form of the velocity distribution 



is a monoscale function of the well pronounced cusp form 
oc exp{ — |w|}. In the case of strong dissipation (ct w 1.0) 
the situation is opposite. The system with low density 
(l — 50), as previously, is characterized by an extremely 
wide velocity distribution, its width is about 10. However, 
now its form deviates substantially from the Gaussian one. 
For the corresponding ensemble with high density {I = 5) 
the velocity distribution is Gaussian with width about 
1. The latter, nevertheless, is much larger then the same 
width in the absence of dynamical traps. 

These features of the velocity distribution character- 
izes the cooperative behavior of particles rather then their 
individual dynamics. In other words, there should be strong 
correlations in the motion of not only neighboring parti- 
cles but also distant ones. Therefore the velocity varia- 
tions responsible for the formation of such distributions 
describe in fact the motion of mutliparticle clusters. To 
justify this we apply to the middle column windows in 
Fig. |31 They visualize some typical fragments of the time 
patterns formed by the velocities of individual particles. 
When the mean particle density is low {I = 50), these pat- 
terns look like a sequence of fragments {wq} inside which 
the particle velocity varies in the vicinity of some level Uq, . 
The values {va} are rather randomly distributed inside a 
certain region of thickness 1/ ~ 10 in the vicinity of u = 0. 
The continuous transitions between these fragments occur 
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Fig. 5. The distribution functions of the particle velocities and the characteristic time patterns formed by the velocity variations 
of 500-th particle. Dynamics of the 1000-particle ensemble with low (/ = 50, label 1) and high {I = 5, label 2) mean density 
and weak (cr « 0.1) and strong (cr « 1.0) dissipation was implemented for the calculation time up to 8000 time units to make 
the obtained distributions stable with respect to time increase. The right windows visualize the time patterns formed by 200 
paths of particle motion during 1000 time units and chosen in the middle of the given ensemble. The other used parameters 
are the noise amplitude e — 0.1, the trap effect measure A = 0.1, the small regularization friction coefficient ctq — 0.01 and the 
regularization spatial scale ho = 0.25. 



via sharp jumps. The typical duration of these fragments 
is about T ~ 100, which enables us to regard them as 
long-lived states because the temporal scales of individ- 
ual particle dynamics are about several units. Moreover, 
these long-lived states can persist if only a group of many 
particles moves as a whole because the characteristic dis- 
tance L individually traveled by a particle involved into 
such state is about L ^ VT 1000 > /. 

The spatial structure of these cooperative states is vi- 
sualized in Fig. 13 right column windows. They depict time 
patterns formed by paths {xi{t)} of 200 particles of du- 
ration about 1000 time units. These particles were chosen 
in the middle part of the 1000-particle ensembles with low 
density. For high density ensembles such patterns also de- 
velop but are not so pronounced. As is seen a large num- 
ber of different mesoscopic states formed in these systems. 
They differ from one another in size, the direction of mo- 
tion, the speed, the life time, ets. Moreover, the life time 
of such a state can be much longer then the character- 
istic time interval during which particles forming it cur- 
rently will belong to this state individually. Besides, the 
found patterns could be classified as hierarchical struc- 
tures. Some relatively small domains formed by coopera- 
tive motion of individual particles in their turn make up 
together larger superstructures. In other words, the ob- 
served long-lived cooperative states have their "own" life 
independent, in some sense, of the individual particle dy- 



namics. The latter properties are the reason for regarding 
them as certain dynamical phases arising in the systems 
under consideration due to the dynamical traps affecting 
the individual particle motion. The term "dynamical" has 
been used to underline that the complex cooperative mo- 
tion of particles is responsible for these long-lived states, 
without the continuous particle motion such states cannot 
exist. 

The obtained results are summarized in the following 
section. 

4 Conclusion 

A rather simple model of an oscillator chain, a one-dimen- 
sional particle ensemble, with dynamical traps and addi- 
tive white noise has been considered. It should be noted 
that the regular "force" governing the individual dynamics 
of particles has no stationary points except for one match- 
ing the system homogeneous state. The latter is locally 
stable for all the possible values of the system parame- 
ters. Nevertheless, as has been demonstrated numerically, 
the sufficiently strong dynamical trap effect accompanied 
with white noise gives rise to a wide variety of anomalous 
and cooperative phenomena. 

In particular, the local symmetry of particle arrange- 
ment (described by variables © called the symmetry dis- 
tortion) can exhibit kinetic phase transitions. Depending 
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on the mean particle density and the dissipation rate the 
distribution function of the symmetry distortion takes ei- 
ther a bimodal form or a twoscale unimodal form, with 
the latter possessing extremely sharp spike (Fig. The 
distortion rate distribution also either is characterized by 
two scales or is of the cusp form oc exp{ — ji^j} smoothed, 
naturally, inside a narrow transition region. 

The cooperative phenomena arising in these system 
have been studied analyzing the velocity distributions and 
visualizing some time pattern formed by the particle dy- 
namics. When the mean particle density is sufficiently 
low the velocity distributions are characterized by the ex- 
tremely large widths. For the system with high density 
and low dissipation the velocity distribution function takes 
also a quasi-cusp form. The visualized time pattern of the 
velocity dynamics of a single particle has demonstrated 
the presence of the long-lived states. They look like a se- 
quence of fragments {wq} within which the particle ve- 
locity varies in the vicinity of some level Va continuously 
joined by sharp jumps in the particle velocity. The velocity 
levels are rather uniformly distributed inside a wide inter- 
val and the life time of these fragments exceeds essentially 
the time scales of the individual dynamics of particles. It 
has been shown that such long-lived states can persist if 
only multiparticle clusters moving as a whole are formed. 

The visualized patterns made up of paths of 200 par- 
ticles have shown also the presence of such cooperative 
structures. Moreover it has become clear that these long- 
lived states persist independently in some sense of the 
individual dynamics of particles forming them currently. 
In other words, the life time of such a state can exceed 
substantially the time interval during which the particles 
forming it at a current time belong to it. Keeping the lat- 
ter in mind we refer to them as to dynamical states. These 
states in turn can form superstructures, so the observed 
patterns are classified as hierarchical structures. 
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